Nematic liquid crystal dynamics under applied electric fields 
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In this paper we investigate the dynamics of hquid crystal textures in a two-dimensional nematic 
under applied electric fields, using numerical simulations performed using a publicly available Liquid 
CRystal Algorithm (LICRA) developed by the authors. We consider both positive and negative 
dielectric anisotropies and two different possibilities for the orientation of the electric field (parallel 
and perpendicular to the two-dimensional lattice). We determine the effect of an applied electric 
field pulse on the evolution of the characteristic length scale and other properties of the liquid crystal 
texture network. In particular, we show that different types of defects are produced after the electric 
field is switched on, depending on the orientation of the electric field and the sign of the dielectric 
anisotropy. 



I. INTRODUCTION 

Defects in the nematic liquid crystalline phase may ap- 
pear when the symmetry of the isotropic phase is broken 
in a isotropic-nematic phase transition, or as a topolog- 
ical response to imposed boundary conditions jT]. Con- 
fined nematics is a stimulating field of recent research 
[2H5] due to the possibility of different technological ap- 
plications. The similarity with defect formation in the 
early universe has also turned liquid crystals into mini 
laboratories for testing the formation and dynamics of 
cosmic topological defects [6-10 . It has also led one of 
us (FM) and co-workers to explore the optical properties 
of defects in a liquid crystal in a geometric setting very 
much like their cosmic analogues |llffH] . 

The complexity of the (dis)organization of nematic 
molecules in the presence of defects limits severely what 
can be analysed using analytical tools. Hence, computer 
simulations are an important complement in the study 
of the evolution of liquid crystalline systems with de- 
fects. Lattice simulations have been used extensively in 
the study of liquid crystals (see, for example, |[T5]-[T2]) 
since the pioneering work of Lebwohl and Lasher j20| . 
An order parameter is usually attached to each point of 
the lattice with orientational degrees of freedom only, so 
that the liquid does not flow but each molecule is free to 
orient itself in three-dimensional space. The fact that a 
relatively large number of particles enter the simulation 
is particularly useful to study collective behavior such 
as phase transitions and systems presenting large scale 
fluctuations and defects. 

From the industry point of view, defects may be quite 
relevant. They are a problem for liquid crystal displays 
since they scatter light, but they are unavoidable in the 
switching of the cells. Consequently, it is interesting, 
both from scientific and technological points of view, to 
study how applied external electric fields infiuence the 
formation and dynamics of defects in nematics. External 
electric or magnetic fields couple to the director adding 



an extra ingredient to an already complex problem |2H - 

In this work we simulate the dynamics of defects in 
a two-dimensional nematic in the presence of an applied 
electric field, after quenching from the isotropic phase. 
The outline of this paper is as follows. In Sec. II we 
present the Landau-de Genes theory which we use to 
model the relaxational dynamics of nematic liquid crys- 
tals under applied electric fields. In Sec. Ill we report on 
a publicly available Liquid CRystal Algorithm (LICRA) 
developed by the authors and describe the specific nu- 
merical implementation used in this paper. In Sec. IV 
the results are presented, considering different possibili- 
ties for the orientation of the electric field and sign of the 
dielectric anisotropy. We conclude in Sec. V. 



II. LANDAU-DE GENNES THEORY 

The orientational order of a nematic liquid crystal is 
described by a symmetric traceless tensor, Q, whose com- 
ponents are given by ^Ij. 

= -l^Sai3+ (fuf{r,U,t)UaUi3, (1) 

where /(r, u, t) is the the molecular orientational distri- 
bution function at the position r and time t, measuring 
the number of molecules with orientation defined by the 
unit vector u with components Uq. The order parameter, 
Q, can be written as a function of the principal values 
S and T and the eigenvectors n (the director), 1 (the 
codirector) and m = n x 1. Its components are given by 
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III. NUMERICAL IMPLEMENTATION 
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XaP = ^ {XaP + Xj3a) — Ti^afl^jj . (3) 

Defining 9 and (p such that cos 9 = un, sin 9 cos </) = u-1 
and sin 9 sin (f> — u ■ m, the principal values S and T can 
be written as 

S = (cos^e)-^, (4) 

T = (sin2 6'cos20). (5) 

The inner product is defined by a ■ h = a-yby where we 
use the Einstein notation for the summation of repeated 
indices. Also, note that in Eq. (hj) (Pu = s\v? 9 dO d(f>. 
The free energy is given by 



F^ j d''v{Fi+Fii + Fui) , 
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where Ti is obtained from a local expansion in powers of 
rotationally invariant powers of the order parameter 

B ^ C ^2 D 

3 ^^ 4 ^2 ^ 6 
Tii is a non-local contribution arising from rotationally 
invariant combinations of gradients of the order parame- 
ter 

J-H = — daQpj daQfjy + — daQal3 djQpj , (8) 

and J^iii accounts for the effect of an external electric 
field E 
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where te = 2Ae and Ae is the dielectric anisotropy. 

Static equilibrium can only be reached for a minimum 
value of the free energy {5F/5Q^i, = 0). However, the 
evolution of the order parameter Q from a given set of 
initial conditions is not fully specified by the free energy 
functionals, and further assumptions have to be made 
on how the minimization process will take place. In the 
absence of thermal fiuctuations and hydrodynamics fiow, 
the time evolution of the order parameter is given by 

6F 
Qa^{r,t)^~T^p,,,——. (10) 

Here the dot represents derivative with respect to the 
physical time, t, and the tensor 
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= thus 



ensuring that the order parameter Q remains symmetric 
and traceless. In the following we shall assume that the 
kinetic coefficient, F, is a constant. 

The equation of motion can be written more explicitly 
as 

Qc^li = -T[{A + CQ2)Qap + iB + DQ3)Q^yQ^p 

— RlQaP,^^ — R2 Ql3'r,ja — Te EaEp ] , (12) 

where a comma denotes a partial derivative. 



We solve Eq. ( 12 I using a standard second-order fi- 



nite difference algorithm for the spatial derivatives and a 
second order Runge-Kutta method for the time integra- 
tion. All the simulations are performed on a two dimen- 
sional 256^ grid (in the xy plane) with periodic bound- 
ary conditions. The order parameter Q has 5 degrees of 
freedom associated with S, T, n, 1 (n accounts for two 
degrees of freedom). In the simulations, the initial condi- 
tions for S and T are randomly generated, at every grid 
point, from uniform distributions in the intervals [0, 1[ 
and [0, S[, respectively. The director, n, is also randomly 
generated, at every grid point, using the spherical vector 
distribuitions routines in the GNU Scientific Library and 
the codirector, 1 , was calculated by randomly choosing 
a direction perpendicular to n. 

We start our simulations with the electric field 
switched off, switching it on at a given time, ion = 400, 
and disconnecting it at a later time, t^s = 600. We 
consider two different possibilities for the orientation of 
the applied electric field, perpendicular or parallel to the 
plane of the simulation. In each case we consider situa- 
tions with either positive (Ae > 0) or negative (Ae < 0) 
dielectric anisotropics. The parameters used in the sim- 
ulations are Ax — Ay = 1, At = 0.1 and A = —0.1, 
B = -0.5, C = 2.67, Ri = 1.0, i?2 = 0.0, T = 1.0, 
corresponding to a uniaxial nematic phase |19j). The 
amplitude of the applied rectangular electric field pulse 
is taken to be ^ = 0.025/ i/JAef. 

In order to determine the evolution of the characteristic 
scale of the network we numerically calculate the function 



F(k,t) 
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after setting the infinite wavelength mode Qap (0) to zero. 
Here Qa/s (k) is the Fourier transform of Qap (r) (k is the 
wavenumber). The characteristic scale, L, can then be 
defined as 



1 
Z2 



= {k')=Y,k'Pik,t)/J2PiKt). 



(14) 



The set of C codes used to perform the simulation are 
available at http://faraday.fc.up.pt/licra along with 
Matlab/Octave routines that are used to generate the 
results visualization. This software is free and we have 
named it LICRA (Liquid CRystal Algorithm). LICRA 
can be used to perform numerical simulation of the dy- 
namics of liquid crystal texture networks in both two- and 
three-dimensional nematics. LICRA also accounts for the 
application of an external electric field on the system. 
Simulation parameters may be modified in the file licra.h 
or directly on the code which is properly commented for 
this effect. The GNU Compiler Collection (gcc) and both 
the Fastest Fourier Transform in the West (FFTW) and 
the GNU Scientific Library (GSL) are required in order 
to compile the source code. 
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FIG. 1: Topological defects in a coarsening nematic liquid 

§stal before the application of an external electric field. Fig. 
represents the value of sin^(2x) at every grid point (note 
that both half-integer and integer disclinations appear in the 
simulations). Fig. lb shows the value of the parameter S 





FIG. 2: Similar to Fig. |la| but now after the application of an 
electric field pulse perpendicular to the lattice plane. The fig- 
ure shows that, in the case of a positive dielectric anisotropy 
(Ae > 0), only integer (±1) disclinations remain in the simu- 
lation, after the application of the external electric field. 



The code outputs contain the results for the order pa- 
rameter S and T, the director and codirector components 
at each grid point and the angle between the projection 
of the director onto the x — y plane and the x axis. The 
visualization routines generate the intensity of the trans- 
mitted light after passing thought a nematic liquid crys- 
tal sample with crossed polarizers, the parameters S and 
T, and the evolution of characteristic length scale with 
time. 



IV. RESULTS 

Fig. 1 represents a texture network in a coarsening 
nematic liquid crystal just before the application of an 
external electric field at i = ton. The value of sin'^(2x) is 
plotted in Fig. |la| where x is the angle of an arbitrary 
axis on the plane of the simulation with the projection 
of the director n onto the same plane. Topological de- 
fects of both half-integer (±1/2) and integer (±1) charge 
appear in the simulations, corresponding to the meet- 
ing points of two and four dark brushes, respectively. 
Schlieren patterns similar to the ones shown in Fig. ( la I 



may be obtained experimentally using a linear polarized 
light microscope with a liquid crystal sample placed be- 
tween crossed polarizers. 

Fig. |lb| shows the value of S at every grid point in- 
side the dashed region as well as the projection, onto the 
two dimensional grid, of the director n. In Fig. [Tc] the 
value of the parameter T is also shown for the same re- 
gion. There is a rapid variation of the parameters 5" and 
T at the half-integer charge defect cores while no signifi- 
cant variation occurs in the case of integer charge defects. 
The non-zero values of T are due to the biaxiality inside 



at every grid point inside the dashed region, as well as the 
projection of the director n onto the two dimensional grid. In 
Fig. [Tclthe value of the parameter T in same region is given. 




(a) 
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(b) 

FIG. 3: Similar to Figs. |la| and [Tb] except that now a negative 
dielectric anisotropy (Ae < 0) was considered (again with an 
electric field pulse perpendicular to the lattice plane). In this 
case, only half-integer (±1/2) disclinations survive, after the 
application of the external electric field. 



half-integer defect cores PSJ. It is also clear that integer 
charge defects are associated with regions where the pro- 
jection of the director n, onto the two-dimensional grid, 
is very small. 

Figs. 2, 3 and 4 represent a texture network in a coars- 
ening nematic liquid crystal after the application of an 
external electric field (i = 800 (> toff))- Fig. l2]shows the 
value of sin (2^) at every grid point after the applica- 
tion of an electric field pulse perpendicular to the lattice 
in the case of a positive dielectric anisotropy (Ae > 0). 
After the field is switched on, the molecules (and the di- 
rector n at each grid point) tend to become aligned par- 
allel to the field, which results in the rapid suppression 
of half-integer disclinations (only integer ones remain in 
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FIG. 4: The figure represents the order parameter S as well 
as the projection of the director, n, onto the two-dimensional 
grid, after the application of an electric field pulse parallel to 
the lattice plane. A negative dielectric anisotropy (Ae < 0) 
was considered in this case. All the defects appearing in the 
simulations correspond to half-integer twisted disclinations. 
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FIG. 5: The evolution of the correlation length, L{t), in four 
different situations: (a) without external electric field, (6) 
with an electric field perpendicular to the lattice plane and 
Ae > 0, (c) the same as (6) but with Ae < and (d) with 
an electric field parallel to the lattice plane and Ae < 0. The 
electric field is connected when t — tan ~ 400 and discon- 
nected when t = ioff — 600. 



the simulations) . Accordingly, S is roughly constant and 
T is very close to zero at all network points. In this case 
the projection of the director n, onto the two-dimensional 
grid, is very small at every grid point. Apart from the 
modification of the background values of S and T there is 
no significant modification to the network evolution when 
the electric field is switched off (this can be confirmed in 
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FIG. 6: Figs. |6a| and |6b| represent the evolution of the pa- 
rameters 5* and T for the case of an homogeneous director, 
n, initiaUy aligned parallel to the lattice plane considering, 
respectively, either a positive (Ae > 0) or negative (Ae < 0) 
dielectric anisotropy. In both cases an electric field perpen- 
dicular to the plane is connected for t = ion = 400 and dis- 
connected for t = tos = 600. 



Fig. [5]). 

Figs. [3a| and [3b| are analogous to Figs. [Ta| and [lb] ex- 
cept that, in this case, a negative dielectric anisotropy 
(Ae < 0) was considered (again with an electric field 
pulse perpendicular to the lattice plane). After the elec- 
tric field is switched on the molecules (and the director 
n) tend to become perpendicular to the electric field, 
which results in the rapid suppression of integer disclina- 
tions (only half-integer ones remain in the simulations). 
As expected there is a rapid variation of the parameter 
S (and T) at the half-integer charge defect cores, which 
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FIG. 7: Figs. |7a| and [Tb] represent the order parameter 5* at 
every grid point at the times corresponding to the maximum 
and minimum of S in Fig. |6b| after the application of an 
electric field perpendicular to the lattice plane, for a negative 
dielectric anisotropy (Ae < 0). 



can be seen in Fig. |3b[ 

If the electric field is applied in a direction parallel 
to the plane of the simulation and Ae > then the di- 
rector n at each grid point will become aligned parallel 
to the electric field and both integer and half-integer de- 
fects are suppressed. The most interesting situation, with 
the electric field parallel to the lattice plane, is that of 
a negative dielectric anisotropy (Ae < 0). In this case, 
the director n at each grid point becomes perpendicular 
to the electric field generating half-integer twist discli- 
nations [131 |5S] . The value of sin^ (2%) is equal to zero 
at every grid point but the defects can be seen in Fig. [4] 
which shows the order parameter S at various grid points 
as well as the projection of the director n onto the lattice 
plane. All the defects appearing in the simulations corre- 



spond to half-integer twisted disclinations which can be 
identified by the values of the projections of n around 
them. 

The evolution of the characteristic length scale, L{t), is 
shown in Fig. [5] in units of the simulation box size. Four 
different situations have been considered: (a) without 
external electric field, (6) with an electric field perpen- 
dicular to the lattice plane and Ae > 0, (c) the same as 
(6) but with Ae < and (d) with an electric field parallel 
to the lattice plane and Ae < 0. In each case an average 
over 100 different simulations with random initial con- 
ditions considering has been made. The electric field is 
connected at the time t — ton — 400 and disconnected at 
t = tos — 600. In all cases the characteristic length scale 
evolves roughly as L^ ex t but the slope is rapidly mod- 
ified when the electric field is switched on (or switched 
off in cases (c) and (d)). 

Figs. [6a| and [6b| show the evolution of the parameters 
S and T for the case of an homogeneous director, n, 
initially aligned parallel to the lattice plane considering, 
respectively, either a positive (Ae > 0) or negative (Ae < 
0) dielectric anisotropy. In both cases an electric field 
perpendicular to the plane is connected for t = ton — 
400 and disconnected for t = toS = 600. As expected. 
Fig. [6a|shows that the system develops a large biaxiality 
right after the electric field is switched on since, initially, 
there are two distinct preferred directions defined by the 
director and the external electric field. After some time, 
the director becomes aligned with the electric field (Ae > 
0), T is again driven towards zero and S tends to larger 
value than before the electric field is switched on (this 
is expected since the dispersion of molecular orientations 
around n, is reduced due to the external electric field). 
After the electric field is switched off at i = ioff , S returns 
to the value it had before ton- 

For Ae < 0, Fig. [6b] shows that again the system 
develops a large biaxiality right after the electric field 
is generated but, in this case, the biaxiality does not go 
away before toff since the two distinct preferred directions 
defined by the director and the external electric field will 
remain perpendicular to each other. Right after the elec- 
tric field is switched on, the value of S increases sharply, 
before moving to a value smaller than before ton- The 
sharp growth of S right after ton is due to an initial de- 
crease of the dispersion of molecular orientations around 
n due to the presence of the external electric field. In 
the subsequent evolution of S, before toff, the molecules 



are forced into the plane by the perpendicular electric 
field. This accounts for the fact that the value of S'(ton) 
is slightly smaller than iS'(toff). Right after toff the value 
of T decreases sharply towards zero, which is responsi- 
ble for corresponding sharp decrease of S. After that, S 
increases again returning to the value it had before ton- 
Figs. [7a| and [7b| represent the parameter, S, at every 
grid point, at the times corresponding to the maximum 
and minimum of S in Fig. |6b[ after the application of an 
electric field pulse perpendicular to the lattice plane, for 
a negative dielectric anisotropy (Ae < 0). They confirm 
the background evolution shown in Fig. |6b| and they also 
reveal a dramatic change on the characteristic scale of 
defect cores which is associated with the behaviour of 
the characteristic scale shown in Fig. [slfor cases (c) and 
(d). 



V. CONCLUSIONS 

In this paper we investigated the effect that an applied 
electric field pulse has on the dynamics of a liquid crystal 
texture network, using numerical simulations performed 
with LICRA, a publicly available Liquid CRystal Algo- 
rithm developed by the authors. We have demonstrated 
that the presence of an electric field can lead to signifi- 
cant changes of the molecular orientational distribution 
function, including the background value of the param- 
eters S" ad T as well as their profile at the defect cores. 
We have further shown that an electric field pulse may 
lead to a defect network whose properties depend on the 
orientation of the applied electric field and the sign of the 
dielectric anisotropy. Depending on the case, topological 
defects with integer or half integer topological charges 
remain after the electric pulse has been applied. In sum- 
mary, our results describe how an applied electric field 
can be used to as control the nature of the defects in a 
liquid crystal texture network. 
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